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1.  PROGRAM  SUMMARY 


The  Space  Science  Group  of  the  Institute  of  ^e^physics  and  Planetary 
Physics  at  the  University  of  California,  Los  Angeles,  is  studying  the 
causes  of  geomagnetic  activity  as  manifested  in  ULF  waves,  magnetospheric 
substorms,  and  magnetic  storms.  Our  research  program  utilizes  a  large  body 
of  magnetic  data  acquired  by  spacecraft  in  the  solar  wind,  magnetosphere, 
and  on  the  earth's  surface,  as  well  as  correlative  information  provided  by 
other  instruments  at  the  same  locations .  i^The  primary  goal  of  our  research 
is  to  understand  the  physical  phenomena  which  transfer  energy  from  the 
solar  wind  to  the  magnetosphere,  which  transport  this  energy  within  the 
magnetosphere,  and  finally  which  lead  to  its  dissipation  in  the  ionosphere 
and  its  return  to  the  solar  wind. 


Our  current  research  projects  include  three  distinct  topics:  solar 
wind'magnetosphere  coupling,  triggering  of  substorm  expansion  onsets,  and 
generation  of  ultra  low  frequency  waves.  Our  investigation  of  solar 
wind-magnetosphere  coupling  is  concerned  with  establi.shing  how  the  rate  of 
energy  transfer  into  the  magnetosphere  depends  on  bulk  parameters  of  the 
solar  wind.  \  Also,  it  involves  a  study  of  how  this  energy  is  dissipated 
internally  b^the  auroral  electro  jets,  the  ring  current,  and  plasmoids  in 
the  magnetotai^.  Our  work  utilizes  linear  prediction  filters  to 
characterize  way  in  which  various  magnetic  indices  are  controlled  by 
the  solar  wind. J 


The  second  topic,  triggering  of  substorm  onsets,  is  concerned  with 
understanding  the  mechanism  whereby  solar  wind  energy  stored  in  the  tail 
lobes  is  released  ioto  the  inner  magnetosphere.  In  the  model  we  favor, 
this  release  is  caused  by  reconnection  of  lobe  magnetic  field  at  a 
near-earth  neutral  line.  Our  present  research  is  concerned  with 
demonstrating  that  this  process  is  highly  localized,  transient,  and  at 
times  involves  multiple  neutral  lines. y  Also,  we  are  trying  to  determine 
whether  this  process  is  triggered  by  ^temal  changes  or  requires  an 
external  stimulus. 


The  third  topic  is  the  study  of  ultra  low  frequency  waves  (ULF) 
(1-1000  mHz) .‘^These  waves  are  generated  by  a  variety  of  mechanisms,  some 
of  which  depend  directly  on  the  solar  wind,  while  others  are  a  consequence 
of  changes  associated  with  substorms  and  storms.  The  subject  of  our 
current  research  is  a  ULF  phenomenon  occurring  on  the  dayside  of  the  earth 
and  referred  to  as  Pc-3  pulsations  (10-4Ss  period).  These  pulsations  are 
apparently  gene-ated  by  fie Id- line  resonances  driven  by  vaves  produced  by 
the  solar  wind  interaction  with  the  magnetosphere.  Our  past  work  shows 
that  these  waves  may  be  used  to  indirectly  measure  the  total  mass  density 
of  the  plasma  in  the  magnetosphere,  enabling  os  to  study  changes  in 
magnetospheric  plasma  associated  with  substorms  and  storms . 


2.1  SOLAR  WIND--MAGNETOSPHERE  COUPLING 


A  major  goal  of  our  current  research  program  is  the  development  of 
methods  of  predicting  geomagnetic  activity  from  measurements  of  the  solar 
wind.  Practically,  such  predictions  would  provide  users  of  geomagnetic 
indices  with  sufficient  warning  that  they  could  modify  their  equipment  or 
procedures  to  avoid  problems  caused  by  severe  substorms  or  major  storms. 
Scientifically,  such  predictions  would  provide  insight  into  the  physical 
processes  which  generate  magnetic  activity.  Recently,  this  area  of 
research  has  come  to  be  known  as  "solar  wind-magnetosphere  coupling."  Our 
investigations  of  this  subject  are  discussed  in  the  following  subsections. 


2.1.1  Oimansional  analysis  of  energy  transfer 

A  major  unsolved  problem  of  solar-terrestrial  physics  is,  determining 
the  functional  dependence  of  the  rate  of  solar  wind  energy  transfer  into 
the  magnetosphere  on  interplanetary  and  magnetospheric  parameters. 

Perreault  and  Akasofu  (1978)  suggested  this  dependence  was  of  the  form 
U.^  »  e(t)  where  U.^  is  the  total  rata  of  magnetospheric  dissipation  and  e  is 
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proportional  to  VB  sin  9/2.  Since  U.^,  depends  both  on  the  rate  of  change 

of  Dst  (ring  current  injection  and  decay)  and  on  AE  (particle  precipitation 
and  Joule  heating)  it  would  be  surprising  if  the  relation  were  so  simple. 

As  discussed  in  Section  2.1.2  below,  our  work  with  Ost  demonstrated  that 
the  rate  of  change  of  Dst  is  quite  predictable  from  the  solar  wind  electric 
field  (Burton  et  al.,  1975).  Also  we  recently  established  that  the  AL 
index  is  more  closely  related  to  the  electric  field  than  it  is  to  s  (Clauer 
et  al.  ,  1981) . 

A  method  for  determining  which  solar  wind  parameter  is  the  best 
predictor  of  a  given  index  of  energy  dissipation  has  been  suggested  by 
Vasyliunas  et  al.  (1982).  This  method  is  based  on  the  principle  of 
dimensional  similitude  which  requires  that  any  magnetospheric  dissipation 
parameter  (P)  be  proportional  to  the  solar  wind  kinetic  energy  flux 
multiplied  by  some  function  of  the  dimensionless  parameters  involved  in  the 
coupling.  For  an  electromagnetic  interaction  (reconnection)  these 
parameters  include  the  orientation  angle  of  the  IMF  and  the  solar  wind  Mach 
number.  Thus, 

P  =  pv"z^jJ(mJ,9) 


where  p  =  density,  V  =  velocity,  =  distance  to  subsolar  point.  A 
simple  approximation  to  F  is  of  the  form 


=  (M^)’®G(9) 


which  leads  to  the  relation 


The  preceding  relation  can  be  written  in  the  form 

log  (P/p^/V^^)  =  a  log  (m2^)  +  log  G(9) 

Providing  the  solar  wind  parameters  are  constant  and  dissipation  has 
reached  equilibrium,  observations  may  be  separated  into  different  data  sets 
corresponding  to  various  IMF  orientations.  Then  a  plot  of  the  L.H.S.  of 
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the  preceding  relation  versxis  M  ^  should  be  a  straight  line  with  slope 
corresponding  to  the  optimum  coupling  exponent. 

We  have  applied  this  method  to  both  the  total  energy  parameter  and 

the  electro jet  index  A£.  Figure  1  illustrates  the  results  when  the  IMF 
orientation  was  selected  in  the  range  -90*’  £  9  ^  -75*’.  A  least  square, 
linear  fit  to  the  data  gives  a  slope  of  0.41  and  a  correlation  coefficient 
of  0.68.  Similar  results  were  obtained  for  all  angles  up  to  ■*‘30*’,  with  a 
slope  0.5  being  most  probable.  For  the  A£  index,  a  slope  of  0.3  is  most 
probable.  These  results  imply  that  the  optimum  coupling  parameter  for  U.p 

is 


U  «  a  (pV^)^^^(VB)G(0) 

This  parameter  is  very  close  to  the  electric  field  parameter  used  by  Burton 
et  al.  (1975)  and  quite  different  from  the  epsilon  parameter  advocated  by 
Akasofu  (1981b)  and  Kan  and  Akasofu  (1982).  For  A£  we  obtain 

AE  .  pl/V/V/2 

This  parameter  is  closer  to  electric  field  than  is  epsilon  but  is  more  like 

the  V^B  dependence  found  empirically  by  several  authors  (e.g.  Murayama, 
1982) . 

A  paper  describing  these  results  was  presented  at  an  international 
symposium  in  Graz,  Austria  (A-94) ,  and  a  manuscript  has  been  published  in 
the  proceedings  (B-69). 


ESTIMATION  OF  COUPLING  EXPONENT 


Figure  1:  Sunmary  of  dimensional  analysis  results  for  UT  parameter.  The 
ordinate  displays  the  optimum  exponent  to  be  used  in  the  solar 
wind-magnetosphere  energy  coupling  equation 


The  abscissa  denotes  different  ranges  of  the  IMF  orientation 


2.1.2  ,  Predicting  the  storm  time  ring  current 

Magnetic  storms  are  disturbances  in  the  earth's  magnetic  field 
produced  by  the  drift  of  particles  in  a  large  doughnut -shaped  volume 
circling  the  earth.  These  particles  are  injected  into  trapped  orbits  from 
the  tail  by  time  varying  magnetospheric  electric  fields,  and  from  the 
ionosphere  by  field  aligned  potential  drops.  They  are  lost  from  these 
orbits  by  precipitation  into  the  atmosphere  through  charge  exchange  and 
wave-particle  interaction.  Because  intense  substorms  accompany  the  growth 
of  this  "ring  current"  it  has  been  cxistomary  to  attribute  particle 
injection  to  these  substorms  (Davis  and  Parthasarathy ,  1967).  .More  recent 
work  by  Perreault  (1974),  and  by  ourselves  (Burton  et  al.,  1975),  however, 
shows  that  the  Ost  index  (a  linear  measure  of  the  energy  in  the  ring 
current)  is  very  closely  related  to  the  solar  wind  electric  field.  As 
discussed  by  Harel  at  al.  (1981),  this  suggests  that  ring  current  injection 
is  more  properly  considered  a  "convection"  phenomenon  driven  by  the  solar 
wind  than  a  "substorm"  phenomenon. 

We  have  attemptad  to  determine  which  of  these  is  more  important  using 
the  technique  of  linear  prediction  filtering.  Applying  the  same  procedures 
as  in  our  AL  study  (B-44)  and  our  ASYM  study  (8-48)  we  have  calculated  a 
prediction  filter  relating  the  solar  wind  electric  field  and  the  Dst  index. 
Because  the  filter  must  represent  both  rapid  injection  and  slow  decay  it  is 
necessary  to  use  a  very  long  filter  and  the  highest  possible  time 
resolution.  In  addition,  because  the  Dst  index  includes  effects  of 
magnetopause  currents,  it  is  necessary  to  correct  for  movement  of  these 
currents  in  response  to  the  solar  wind  dynamic  pressure.  This  was  done  by 
calculating  a  prediction  filter  relating  Dst  and  dynamic  pressure  and 
subtracting  its  prediction  from  the  observed  index.  This  residual  should 
contain  only  the  effects  of  solar  wind  electric  field. 

Our  results  show  that  the  dynamic  pressure  filter  is  approximately  a 
delta  function,  while  the  electric  field  filter  is  a  pulse  which  rises  to  a 
peak  after  only  one  hour,  and  which  decays  slowly  to  background  sometime 
after  16  hours.  Separately,  the  pressure  filter  accounts  for  12%  of  the 
Dst  variance,  and  the  electric  field  filter  45%,  while  together,  they 
account  for  72%  of  the  variance.  An  example  illustrating  the  quality  of 
the  predictions  is  presented  in  Figure  2.  As  apparent  from  the  comparison 
with  observations  in  the  bottom  panel,  there  are  significant  residuals. 
These  appear  to  arise  from  a  variety  of  sources  including  errors  in 
measurement  of  the  upstream  parameters,  the  definition  of  the  quiet  dav 
used  in  calculating  the  observed  Dst  index,  and  deviations  of  the  decay 
rate  from  the  average  represented  by  the  filter. 

Initial  results  of  this  work  are  contained  in  a  paper  presented  at  the 
COSPAR  meeting  in  Graz  (A-97;  B-68).  A  manuscript  describing  the  complete 
investigation  is  in  preparation. 


2.1.3  Predicting  the  Polar  Cap  Potential 


Our  past  efforts  to  predict  magnetic  activity  have  made  use  of 
magnetic  indices  because  they  are  so  readily  available.  A  potential 
objection  to  this  use  is  that  the  indices  are  only  indirectly  related  to 
the  underlying  physical  processes.  An  opportunity  to  use  more  physical 
parameters  has  been  provided  by  several  recent  workshops  (Kamide  et  al . , 
1982;  McPherron  and  Manka,  1985).  In  these  workshops  a  large  collection  of 
ground  magnetic  observations  have  been  analyzed  to  produce  a  model  of 
"true"  ionospheric  currents,  field  aligned  currents,  ionospheric  electric 
field  and  Joule  heating. 

In  our  study  of  these  data,  we  began  by  calculating  and  comparing 
prediction  filters  for  AU  and  AL  and  their  corresponding  physical 
parameters,  eastward  and  westward  current  strength.  Within  the  constraints 
of  the  limited  amount  of  data  available  for  the  study,  we  cannot 
distinguish  between  the  corresponding  pairs  of  filters.  This  suggests  that 
our  previous  work  based  on  the  indices  alone  has  not  led  to  inappropriate 
conclusions  as  a  result  of  the  imperfect  way  in  which  the  indices  monitor 
the  underlying  physical  quantities. 

Our  results  for  the  pol^r  cap  potential  are  more  interesting,  however, 
because  this  potential  is  used  extensively  by  the  Rica  group  as  input  in  a 
numerical  simulation  of  magnetospheric  convection  (cf . ,  Harel  et  al., 

1981).  In  the  past,  polar  cap  potential  has  been  calculated  from  a  simple 
time* independent  model  of  the  relation  between  polar  cap  potential,  and 
solar  wind  (Reiff  et  al.,  1981).  An  empirical  model  which  takes  into 
account  the  time  delays  in  the  magnetospheric  response  should  be  helpful  in 
improving  these  simulations.  Also,  these  delays  may  provide  insight  into 
the  physical  processes  involved  in  the  interaction. 

Our  results  are  shown  in  Figure  3.  Two  filters  are  required  to 
predict  this  potential.  The  first  is  for  the  rectified  solar  wind  electric 
field,  while  the  second  is  for  the  square  of  the  solar  wind  velocity.  The 
latter,  which  has  the  form  of  a  delta  function,  is  required  to  predict  the 
observed  potential  drop  during  times  of  northward  IMF.  The  form  of  filter 
which  transforms  solar  wind  electric  field  to  polar  cap  potential  is 
difficult  to  perceive  because  of  noise  in  the  filter  estimates.  It  appears 
to  rise  to  a  maximum  at  40  minutes  and  decay  to  zero  after  90  minutes.  The 
steady  state  response  to  a  1  mV/m  electric  field  is  20  keV.  For 
comparison,  the  resting  potential  for  an  extremely  slow  solar  wind 
(v  =  200  km/sec)  and  northward  IMF  would  be  10  keV.  For  more  _ypical 
conditions  (5  mV/m  and  400  km/sec)  the  values  are  100  jteV  and  40  keV, 
respectively.  While  the  effect  of  the  electric  field  is  strongest,  the 
effect  of  the  velocity  (viscous  interaction)  is  very  important. 
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Figure  3:  Filters  which  predict  polar  cap  potential.  Top  panel  shows 
solar  wind  electric  field  filter.  Bottom  panel  shows 
velocitv  filter. 


1 


2.1.4  Application  of  Activity  Predictors  to  CDAW-6  Observations 


Our  empirical  technique  for  predicting  geomagnetic  activity  utilizes 
past  observations  to  define  a  system  response  function  which  is  then 
convolved  with  a  later  input  to  "predict"  the  output.  If  this  input  is  the 
same  as  used  to  define  the  response  function,  the  output  is  not  truly  a 
prediction.  Unless  the  underlying  system  is  approximately  linear  and  time 
invariant,  this  function  is  only  a  "representation"  of  a  particular  set  of 
data.  If,  however,  the  result  of  convolution  with  an  arbitrary  input  is 
reasonably  close  to  the  observed  output,  then  the  response  function  may  be 
considered  a  predictor  as  well. 

In  order  to  test  the  predictive  value  of  the  response  functions 
calculated  in  our  previous  studies,  we  have  performed  "predictions"  for 
two- three  day  intervals  which  encompass  the  events  chosen  for  study  in  the 
Coordinated  Data  Analysis  Workshop,  CDAW-6  (see  Section  2.2  below).  These 
intervals  were  chosen  because  the  CDAW-6  workshops  have  provided  a 
generally  accepted  morphology  with  which  the  residuals  in  our  predictions 
can  be  compared.  Also,  high  time  resolution  magnetic  indices  based  on  a 
large  number  of  stations  were  available  for  comparison  with  the 
predictions . 

The  test  was  carried  out  in  two  stages.  First,  the  52 -day  data  set 
from  1967-68  was  averaged  to  5 -minute  resolution  to  correspond  to  the 
resolution  of  the  CDAW-6  solar  wind  data.  Impulse  response  functions  for 
all  four  indices  were  then  calculated  at  this  resolution  and  plotted. 

Next,  these  functions  were  convolved  with  the  CDAW-6  solar  wind  input  to 
predict  each  index.  Finally,  the  observed  indices  were  plotted  on  cop  of 
the  predictions  for  comparison.  Figure  <*  is  an  example  of  such  a 
comparison . 

In  summary,  we  find  the  AU  index  is  least  predictable  (-20%),  AL  and 
ASYIl  moderately  predictable  (-40%)  and  Dst  most  predictable  (-72%) .  We 
also  find  there  is  virtually  no  difference  between  the  residuals  in  the 
data  used  to  define  the  response  functions  and  the  CDAW-6  residuals 
suggesting  the  response  functions  truly  represent  physical  properties  of 
the  system.  We  attribute  the  residuals  in  the  Dst  index  to  errors  in  the 
input  data,  to  problems  defining  quiet  base  lines  and  finally  to  an 
inability  to  accurately  define  the  decay  rate  of  the  ring  current.  The 
latter  is  probably  a  consequence  of  ionic  composition  depending  on  the 
level  of  activity.  We  attribute  the  residuals  in  AL  to  a  ) ack  of  a 
deterministic  relation  between  the  solar  wind  and  the  exact  time  and  amount 
of  energy  released  from  the  tail  during  a  substorm  expansion.  Our 
inability  to  predict  the  AU  index  is  difficult  to  understand  and  suggests 
the  eastward  eleccrojec,  which  it  measures,  depends  strongly  on  parameters 
other  chan  the  solar  wind  electric  field. 

These  results  have  been  presented  at  a  recent  meeting  in  Graz,  Austria 
(A-97)  and  have  been  published  in  the  proceedings  of  that  meeting  (3-68). 


2.1.5  Roi«s  of  driven  and  loading/unloading  processes  in  substorms 


For  over  a  decade  it  has  been  known  that  the  magnetotail  plays  a  major 
role  in  magnetospheric  substorms.  As  outlined  by  McPherron  et  al.  (1973) 
and  Russell  and  McPherron  (1973),  this  role  is  to  accumulate  energy  from 
the  solar  wind-magnetosphere  interaction  and  to  later  release  this  energy 
in  bursts  called  substorms.  Recently,  an  alternative  model  has  been 
suggested  by  Akasofu  (1981a)  in  which  substorm  phenomena  are  explained  as 
manifestations  of  field  aligned  and  ionospheric  currents  directly  driven  by 
the  solar  wind-magnetosphere  interaction.  In  this  model  the  only  role  of 
the  magnetotail  is  to  determine  the  size  of  the  polar  cap,  and  hence,  the 
location  of  disturbances  in  the  ionosphere.  Akasofu  (1979)  characterizes 
these  alternate  models  as  "loading/ unloading"  systems  and  "driven"  systems 
respectively.  In  numerous  publications  (e.g.  Akasofu,  1981b)  he  has  argued 
that  the  magnetosphere  is  primarily  a  driven  system  and  that  the 
magnetotail  plays  a  passive  role. 

The  importance  of  driven  processes  in  magnetospheric  dynamics  is 
implicit  in  the  loading/ unloading  model  of  substorms.  The  process  of 
energy  storage  involves  merging  of  interplanetary  magnetic  field  with  the 
earth's  field  and  transport  of  dayside  magnetic  flux  to  the  tail  lobes. 
Subsequently,  the  field  lines  are  reconnected  and  returned  to  the  dayside 
completing  a  closed  circuit  of  magnetospheric  convection.  Field  aligned 
currents  generated  by  this  convection  close  through  the  ionosphere  and 
produce  significant  and  observable  effects.  In  quasi-steady  state  there 
processes  are,  in  fact,  the  driven  system.  In  the  loading/ unloading  model, 
however,  it  is  argued  chat  these  processes  bring  about  changes  in  the 
magnetospheric  configuration  which  ultimately  lead  to  the  rapid  release  of 
energy  accumulated  in  the  tail  lobes.  Much  of  this  energy  is  dissipated  in 
the  atmosphere  as  is  the  case  for  the  driven  processes.  The  question  of 
whether  the  magnetosphere  is  primarily  a  driven  or  an  unloading  system  is  a 
question  of  the  relative  magnitude  of  the  energy  dissipation  associated 
with  the  two  types  of  processes.  An  attempt  to  pose  this  question,  and  to 
define  driven  and  unloading  processes  was  made  in  a  manuscript  submitted 
for  publication  last  November  (6-50).  Also,  we  have  attempted  to  answer 
the  question  by  a  careful  study  of  several  substorms  occurring  during  the 
CDAW-6  intervals  (see  subsection  2.2.2  below).  A  manuscript  (B-71) 
describing  this  CDAW-6  study  has  been  accepted  for  publication. 


2.1.6  Substorm  Triggering  by  the  Interplanetary  Magnetic  Field 

In  the  loading/ unloading  model  of  substorms  energy  extracted  from  the 
solar  wind  is  stored  in  the  tail  lobes  and  later  released.  A  variety  of 
mechanisms  for  initiating  this  release  have  been  suggested.  Among  these 
are  solar  wind  pressure  discontinuities  (Kokobun  et  al.,  1977)  plasma  sheet 
thinning  (McPherron  et  al.,  1973),  field  aligned  current  instability 
(Akasofu,  1981a),  and  northward  turnings  of  the  IMF  (Caan  et  al.,  1977). 

The  last  of  these  possibilities  has  attracted  recent  attention  and  Rostoker 
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(1983)  and  Rostoker  et  al.  (1983)  have  published  several  examples  of 
substorms  which  appear  to  have  been  unambiguously  triggered  by  northward 
turnings  of  the  IMF.  These  results  have  led  Rostoker  to  conclude  that  all 
major  substorm  expansions  are  so  triggered  (Rostoker  et  al.,  1983). 

While  it  is  our  own  work  which  led  to  the  identification  of  IMF 
triggering,  (e.g.,  Caan  et  al.,  1977,  1978)  we  do  not  believe  this  is  the 
only  cause  of  onset,  as  it  does  not  explain  why  the  typical  delay  time 
between  a  southward  turning  and  expansion  onset  is  about  60  minutes.  This 
could  happen  if  60  minutes  were  the  typical  duration  of  southward  IMF,  but 
this  is  not  observed.  We  feel  it  is  far  more  likely  that  this  delay 
represents  a  time  scale  inherent  to  the  magnetosphere,  that  is,  the  time 
for  the  tail  current  to  reach  an  instability  threshold. 

We  have  studied  the  question  of  IMF  triggering  by  attempting  to  find 
unambiguous  examples  of  both  IMF  and  internally  triggered  onsets.  In  this 
investigation  we  find  that  there  are  few  such  examples.  There  are  several 
reasons  for  this.  First,  the  IMF  is  seldom  steady  in  direction  and  often 
fluctuates  widely.  At  such  times  an  apparent  correlation  can  always  be 
made.  This  is  particularly  true  because  it  is  not  obvious  when  an  observed 
discontinuity  should  arrive  at  the  earth.  Our  work  shows  that  several 
published  correlations  are  probably  fortuitous  since  careful  analysis  of 
the  propagation  of  the  IMF  discontinuity  to  the  earth  gives  an  arrival  time 
different  from  the  expansion  onset.  Also,  we  find  the  concept  of 
"discontinuity"  ambiguous  as  well.  In  several  examples  of  expansion  onset 
during  what  we  consider  "steady"  IMF,  Rostoker  (personal  communication)  has 
pointed  out  exceedingly  small  changes  in  slope  (i.e.  constant,  to  1  nT/hour 
northward  turning  as  the  probable  trigger.  Figure  5  depicts  such  a 
situation. 

In  view  of  the  foregoing  we  conclude  that  it  is  not  possible  to  decide 
whether  or  not  triggering  by  the  IMF  occurs  with  the  data  available  to  us 
for  this  study.  Our  preliminary  results  were  presented  at  the  December 
1983  AGU  meeting  (A-95)  and  a  draft  manuscript  has  been  completed. 
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2.2  STUDIES  OF  THE  NEAR-EARTH  PLASMA  SHEET 
2.2.1  Dynamics  of  an  isoiatad  substorm 

In  the  past  decade  there  have  been  a  number  of  papers  describing 
detailed  time  histories  of  magnetospheric  substorms  Ce.g.  McPherron  et  al., 
1973;  Nishida  and  Kamide,  1983).  The  primary  goal  of  these  studies  has 
been  to  determine  the  physical  processes  by  which  energy  and  particles  are 
transferred  from  the  solar  wind  to  the  magnetosphere  and  transported  within 
to  final  deposition  in  the  ionosphere.  With  the  introduction  of  the 
Coordinated  Data  Analysis  Workshops  (CBAW)  (Vette  et  al.,  1982),  such 
studies  have  become  increasingly  detailed.  In  CDAW  workshops  data  from 
many  spacecraft  and  ground  observatories  are  combined  in  a  digital  data 
base  and  jointly  studied  by  a  number  of  investigators  from  many  countries. 

Several  years  ago  we  initiated  the  investigation  of  two  intervals  of 
substorm  activity  in  a  sequence  of  workshops  which  have  been  named  CQAW-6. 
In  this  investigation  the  primary  emphasis  has  been  an  isolated  substorm 
which  followed  an  interval  of  magnetic  quiet,  and  during  which  the  ISEE  1/2 
spacecraft  were  inbound  in  the  near-earth  plasma  sheet.  Fortunate 
circumstances  placed  the  spacecraft  on  the  central  meridian  of  the 
disturbance,  and  also  in  the  plasma  near  the  center  of  the  neutral  sheet  at 
the  time  of  expansion  onset. 

Our  major  contribution  to  this  study  has  been  an  overview  of  the  event 
establishing  the  time  history  of  the  substorm,  and  the  general 
correspondence  of  the  observations  to  the  "neutral -line"  model  of 
substorms.  In  our  opinion,  the  gross  behavior  of  the  magnetic  field  and 
plasma  supports  the  idea  that  magnetic  merging  on  the  dayside  enables  the 
solar  wind  to  transport  dayside  magnetic  flux  to  the  lobes  of  the  tail.  A 
concurrent  return  of  flux  exhausts  the  plasma  sheet  as  simultaneously  the 
lobe  flux  increases.  Eventually,  the  near-earth  plasma  sheet  thins  to  the 
point  that  reconnection  begins,  diverting  a  portion  of  the  enhanced  tail 
current  to  the  ionosphere.  When  the  merging  process  encounters  open  field 
lines  of  the  lobe,  a  plasmoid  is  formed  and  pulled  from  the  tail  by 
reconnected  field  lines  of  the  solar  wind,  .'lerging  continues  until 
quenched  by  exhaustion  of  available  energy,  or  by  a  relaxation  of  solar 
wind  stress  (northward  turning  of  IMF) .  At  this  point  the  neutral  line 
moves  rapidly  tailward  expanding  the  plasma  sheet. 

Although  the  gross  behavior  of  plasma  and  field  seems  consistent  with 
this  model,  opponents  of  this  view  have  raised  a  number  of  questions, 
particularly  with  behavior  of  more  energetic  particles  in  the  plasma  sheet 
boiindary  layers.  In  our  opinion,  these  can  be  explained  by  the  formation 
of  multiple  neutral  lines  in  different  azimuthal  and  radial  sectors  at 
different  times  in  the  substorm.  Unfortunately,  these  assumptions  cannot 
be  verified  with  the  available  data,  and  opponents  of  the  neutral  line 
model  remain  unconvinced. 
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Our  description  of  the  dynamics  of  this  substorm  has  been  presented  in 
papers  at  several  meetings  (A-79,  A-85).  A  manuscript  describing  the 
results  has  been  accepted  for  publication  (B-64) .  Our  paper  will  appear  as 
the  lead  paper  in  a  set  of  approximately  20  papers  in  the  Journal  of 
Geophysical  Research.  February  1985. 


2.2.2  Magnmtotaii  anargy  storaga  and  raiaasa  during  substorms 

One  of  the  outstanding  questions  concerning  substorms  is  the  relative 
importance  of  driven  and  loading/unloading  processes.  As  discussed  above 
in  Section  2.1.5,  it  is  now  evident  that  both  processes  contribute  to 
substorms,  however,  which  process  produces  the  most  energy  dissipation 
cannot  be  determined  from  studies  of  magnetic  indices.  According  to 
Akasofu  (1981b)  the  similarity  between  the  solar  wind  energy  parameter  e 
and  the  A£  and  indices  implies  that  the  magnetospheric  system  is 

predominantly  driven.  However,  as  we  have  pointed  out,  if  the  release  of 
energy  from  the  magnetotail  is  a  consequence  of  solar  wind  induced  changes 
in  the  tail,  the  effects  of  the  release  will  also  be  correlated  with  the 
solar  wind. 

One  way  of  addressing  the  question  of  relative  importance  is  to 
compare  physical  measures  of  energy  dissipation  at  times  when  only  one  or 
the  other  of  the  two  processes  is  operative.  An  opportunity  to  make  this 
comparison  was  provided  by  the  CDAW'd  data  base  described  in  section  2.2.1. 
The  large  collection  of  ground  magnetic  data  (Kamide  et  al.,  1983)  made  it 
possible  to  model  ionospheric  currents  and  the  Joule  heating  they  caused. 
Also,  an  accurate  Dst  provided  a  measure  of  dissipation  due  to  ring  current 
injection.  The  same  data  also  provided  the  information  needed  to  define 
substorm  phases,  in  particular  expansion  onsets  which  are  generally  thought 
to  be  caused  by  release  of  energy  from  the  tail.  Finally,  the  availability 
of  satellite  data  made  it  possible  to  determmine  the  state  of  the 
magnetotail . 

Our  study  shows  that  the  time  history  of  Joule  heating  during  a 
substorm  has  two  components .  The  first  begins  soon  after  a  southward 
turning  of  the  IMF,  the  second  shortly  after  the  onset  of  an  expansion 
phase.  In  cases  where  the  IMF  remains  rather  steadily  southward,  the  Joule 
heating  curve  returns  to  the  pre-expansion  level  after  about  one  hour. 

This  behavior  suggests  the  first  component  is  the  driven  component  while 
the  latter  is  the  unloading  component.  Figure  6  illustrates  this  behavior 
for  a  short-  and  a  long-duration  event. 

A  comparison  with  other  data  indicates  that  the  second  component 
occurs  after  expansion  phase  onset,  and  in  coincidence  with  dipolarization 
of  the  field  at  synchronous  orbit  and  a  decrease  in  tail  lobe  field 
intensity. 

From  the  figure,  it  is  evident  that  the  magnitude  of  the  Joule  heating 
rate  after  expansion  onset  is  4-3  times  as  large  as  that  before  or  after. 
Extrapolating  the  pre-onset  values  to  later  times  and  integrating  the  area 
under  this  curve  gives  a  total  heating  due  to  the  "driven”  component. 
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Subtracting  this  from  the  total  integral  gives  the  component  due  to 
unloading.  We  find  that  the  unloading  component  dissipates  2-5  times  as 
much  energy  as  the  driven  component. 


2.2.3  Changes  in  the  plasma  shaat  magnetic  field  during  substorms 

As  part  of  the  CDAW-6  investigation  described  in  Section  2.2.1,  we 
have  carried  out  a  detailed  study  of  the  plasma  sheet  magnetic  field.  The 
motivation  for  this  study  was  the  observation  that  shortly  after  onset  the 
field  at  the  neutral  sheet  was  directed  southward  with  magnitude  ~60  nT  and 
that  it  was  accompanied  by  a  rapid  tailward  flow  of  plasma.  Such  behavior 
is  expected  tailward  of  a  neutral  line  which  in  this  case  must  have  formed 
earthward  of  -14  Re. 

Our  ability  to  observe  this  event  is  a  consequence  of  special 
circumstances  of  dipole,  spacecraft  and  solar  wind  geometry.  At  the  time 
of  the  event,  1055  UT  March  22,  1979,  the  dipole  axis  was  almost  exactly 
perpendicular  to  the  earth-sun  line  so  that  the  neutral  sheet  should  be 
located  in  the  GSM  equatorial  plane,  with  the  spacecraft  ~1  Re  south  of  it. 
Fortunately,  however,  the  southward  turning  of  the  IMF  which  generated  the 
substorm  was  accompanied  by  a  southward  turning  of  the  solar  wind  velocity 
vector.  Later,  the' northward  turning  of  the  IMF  was  accompanied  by  a 
corresponding  northward  change  in  velocity.  These  changes,  combined  with 
waves  on  the  neutral  sheet,  caused  multiple  crossings  of  the  neutral  sheet 
which  were  observed  by  both  ISEE  spacecraft.  Using  time  delays  between  the 
two  spacecraft  we  have  been  able  to  estimate  the  velocity  and  thickness  of 
the  plasma  sheet,  as  well  as  transform  the  magnetic  data  to  an  appropriate 
"tilted"  neutral  sheet  coordinate  system. 

Our  results  show  that  at  the  beginning  of  the  substorm  the  current 
sheet  was  quite  thick.  However,  as  the  growth  phase  progressed  it  rapidly 
thinned.  Immediately  after  expansion  onset  it  became  so  thin  (SIOOO  km) 
that  the  two  spacecraft  were  on  opposite  sides  of  the  sheet  for  almost  15 
minutes.  Furthermore,  the  weakly  southward  IMF  seen  throughout  this 
interval  by  both  spacecraft  was  truly  a  southward  field  and  not  an  effect 
of  tilting  of  the  neutral  sheet.  Finally,  timing  of  the  intense  southward 
field  at  the  two  spacecraft  suggests  that  a  plasmoid  was  ejected  from  the 
near-earth  region  with  a  velocity  comparable  to  the  plasma  streaming 
velocity  (500  km/sec) . 

A  paper  describing  these  observations  has  been  presented  at  two 
meetings  (A-80,  A-86)  and  a  manuscript  is  under  preparation.  Figure  7 
summarizes  our  interpretation  of  changes  in  field  configuration  at  v  .rious 
times . 
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Figure  7;  Interpretation  of  magnetic  field  observations  during  1054  UT 
.March  22,  1979  substora.  Top  left  panel  shows  possible  quiet 
time  plasma  sheet  configuration.  Bottom  left  panel  shows 
configuration  immediately  before  onset.  Upper  right  panel  shows 
configuration  shortly  after  onset. 
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2.2.4  Changes  in  th«  Synchronous  magnetic  field  during  substorms. 

During  the  CDAW-6  substoros  of  March  22,  1979,  two  synchronous 
spacecraft  carrying  magnetometers  passed  through  the  midnight  sector.  As 
we  have  shown  previously  (Kokobun  and  McPherron,  1981)  the  synchronous 
magnetic  field  exhibits  a  "two  phase"  behavior  during  substorms.  During 
the  growth  phase,  the  field  orientation  becomes  increasingly  tail-like, 
while  during  the  expansion  phase  it  becomes  more  dipolar..  Baker  et  al. 
(1978)  have  reported  similar  behavior  from  energetic  particle  anisotropies 
at  synchronous  orbit  and  suggest  that  such  measurements  can  be  used  to 
predict  the  occurrence  of  substorm  expansions.  Nagai  (1982)  has  extended 
these  observations  using  multiple  synchronous  spacecraft  magnetometers  to 
determine  the  spatial  extent  of  the  "tail- like"  and  "dipole- like"  field 
changes .  He  concludes  that  the  former  is  a  result  of  inward  movement  and 
growth  of  the  tail  current,  while  the  latter  is  caused  by  the  rapid 
eastward  and  westward  expansion  of  a  field  aligned  current  wedge  formed  at 
expansion  onset. 

Using  data  from  the  CDAW-6  data  base,  we  have  examined  the  behavior  of 
the  synchronous  magnetic  field  during  the  March  22,  1979  substorms.  We 
find  that  the  earliest  indication  of  expansion  onset  observed  anywhere  in 
the  magnetosphere  was  at  a  synchronous  spacecraft  located  on  a  field  line 
very  close  to  the  meridian  and  latitude  at  which  the  westward  traveling 
surge  formed.  Seventeen  minutes  later,  the  region  of  disturbance  expanded 
and  a  synchronous  magnetometer  30^  east  of  the  first  recorded  a 
dipolarization  of  the  field. 

These  results  support  the  idea  that  the  region  in  which  the  substorm 
expansion  begins  is  close  to  the  earth  and  involves  a  diversion  of  tail 
current  into  the  ionosphere.  Such  a  diversion  is  required  for  the 
formation  of  a  neutral  line,  which  is  the  interpretation  we  made  of  the 
near-earth  plasma  sheet  measurements  (c.f.  Section  2.2.3).  The  results 
also  support  the  suggestion  of  Nagai  that  the  region  of  diversion,  i.e.  the 
neutral  line,  expands  azimuthally  during  the  substorm  expansion  phase.  A 
schematic  illustration  of  this  current  system  is  presented  in  Figure  3. 

A  manuscript  describing  these  results  has  been  submitted  for 
publication  as  part  of  the  CDAW-6  collection  (B-70)  . 


Figure  3;  Schematic  of  currents  causing  expansion  phase  magnetic 

perturbations.  The  east  and  west  edges  expand  away  from  the 
center  as  time  progresses . 


2.3  ULTRA  LOW  FREQUENCY  WAVES 


Ultra  low  frequency  (ULF)  waves  are  observed  throughout  the 
magnetosphere  virtually  all  the  time.  They  are  generated  by  a  variety  of 
mechanisms  such  as  field  line  resonance  and  plasma  instabilities.  Because 
the  wave  velocities,  frequencies,  and  polarization  depend  on  the  magnetic 
field  and  plasma  conditions  in  the  region  of  generation,  they  provide  a 
means  of  indirect  measurement  of  such  quantities.  During  the  past  year  we 
have  continued  our  study  of  these  waves,  emphasizing  Pc-3  pulsations  which 
appear  to  be  generated  in  the  solar  wind. 


2.3.1  Pc-3  magnetic  pulsations 

Pc-3  (pulsations  continuous  15  <  T  <  45  s)  are  a  daytime  phenomenon 
observed  both  at  the  earth's  surface  and  in  space.  Their  occurrence  is 
closely  controlled  by  the  solar  wind  velocity  and  magnetic  field 
orientation  (B-39).  At  synchronous  orbit  their  complex  waveform  can  be 
resolved  into  a  sequence  of  harmonic  frequency  bands,  the  separation  of 
which  gives  the  total  mass  density  distributed  along  a  field  line  passing 
through  the  spacecraft  (B-47).  Observations  with  multiple  spacecraft 
demonstrate  that  every  point  in  space  is  oscillating  with  a  different 
frequency  spectrum  dependent  on  local  conditions  CB-33).  The  magnetic 
polarization  of  these  oscillations  is  always  azimuthal,  never  radial  as 
previously  reported  (B-55),  and  corresponds  to  predictions  made  by  theories 
of  driven  field  line  resonances  (Southwood,  1974).  Their  relation  to  the 
solar  wind  suggests  they  are  generated  outside  the  magnetopause  and 
penetrate  through  the  boundary  to  drive  resonant  field  line  oscillations. 

To  stimulate  the  broad  range  of  harmonics  observed,  the  input  .spectrum  must 
be  very  broad.  A  suggestion  chat  the  period  of  Pc-3  pulsations  is 
controlled  by  the  magnitude  of  the  upstream  solar  wind  magnetic  field 
(Gul'yel'mi,  1973)  does  not  seem  to  be  true  at  synchronous  orbit.  The  only 
observable  IMF  effects  is  a  slight  increase  of  the  amplitude  of  the  highest 
harmonics  with  increasing  field  magnitude  (B-54) . 

During  the  last  year  our  main  activity  in  the  study  of  Pc-3  pulsations 
has  been  to  complete  a  PhD  dissertation  on  this  topic  (Takahashi,  1983), 
and  to  publish  the  papers  summarized  in  the  preceding  paragraph.  .New 
results  contained  in  the  dissertation  and  not  yet  published  are  concerned 
with  the  relative  amplitude  and  frequency  of  harmonics,  and  the  sensitivity 
of  field  line  eigenf requencies  to  boundary  conditions  field  models  and 
plasma  models.  The  observed  distribution  of  the  ratio  of  third  to  second 
harmonic  frequency  is  shown  in  panel  a  of  Figure  9.  The  peak  of  the 
histogram  at  1.56  implies  chat  the  most  probable  distribution  of  plasma 
mass  along  a  field  line  is  proportional  to  the  fourth  power  of  radial 
distance.  Using  this  model  we  have  calculated  numerically  the  expected 
ratio  of  harmonic  amplitudes  as  a  function  of  latitude,  (panel  b).  For 
equal  e.xcitation,  Che  expected  inequality  at  10®  magnetic  latitude  is 

A.  >  A,  >  A.  >  A.  »  A. . 
t  3  6  2  >+ 


Panel  a  of  Figure  10  contains  the  experimental  results  where  it  is  apparent 
that  the  inequality  for  median  amplitudes  is  more  like  >>  A.  -  A^. 

The  major  discrepancy  is  that  A.  and  A^  are  very  much  smaller  that 

predicted.  A  likely  explanation  for  this  is  that  the  input  spectrum 
decreases  at  frequencies  above  the  third  harmonic.  In  fact,  as  noted 
above,  there  is  a  weak  dependence  of  these  harmonic  amplitudes  on  the 
strength  of  the  IMF. 

Another  factor  which  may  contribute  to  this  discrepancy  is  the 
sensitivity  of  the  location  of  nodes  of  the  higher  harmonics  to  the 
distribution  of  plasma  along  field  lines.  For  example,  as  demonstrated  in 
panel  b  of  Figure  10,  the  location  of  the  fourth  harmonic  node  changes  by 
over  10®  depending  on  the  exponent  used  to  characterize  the  plasma 
distribution.  Also,  as  shown  by  the  dashed  line,  the  actual  latitude  of 
the  spacecraft  is  a  fiinction  of  time,  and  hence  it  must  move  through  nodes 
of  higher  harmonics  throughout  the  day. 

The  frequency  analysis  leads  us  to  the  conclusion  chat  our  current 
knowledge  of  the  spacecraft  position  in  the  actual  magnetic  field,  and  the 
mass  distribution  along  field  lines,  are  not  well  enough  known  to  predict 
relative  harmonic  amplitudes. 


Distribution  of  1^/(2 
ATS  6  Azimuthal  Component 
August  1974,  August  1975 


jr  WKB  Approximation 
:6  5  4  3  2  10 

)  I  11  111 

t  Median 


•  • 
•  • 


•  •• 


1.55  1.60 


Q  Oh- 


'x>*Jty-AUGUST.  STS  '  S74-.*aY  STS 


361—  / 


0  2 


4  o  3  O  '2  14 

"4AGNeTIC  UkTITUOe  (OEGSEE) 


izure  9;  Expected  harmonic  amplitudes  at  svnchronous  orbit.  (a.)  Observed 
ratio  of  third  to  second  harmonic  freauencies  are  compared  to 
expectations  for  different  power  law  distributions  of  plasma  mass 
density.  (b)  Calculated  wave  amplitude  as  a  function  of  latitude 
for  different  harmonics  using  most  probable  mass  distribution 
(m  *  4)  . 
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2.4  SOFTWARE  DEVELOPMENT 
2.4.1  Flatfila  software 


For  several  years,  tae  Space  Science  Group  of  the  Institute  of 
Geophysics  and  Planetary  Physics  at  UCLA  has  been  developing  software  to 
support  its  activities  in  management  and  analysis  of  time  series  data. 
During  the  past  year,  we  have  concentrated  on  extending  this  system  from 
the  HP  minicomputer  to  an  APPLE  and  an  IBM  personal  computer.  The  primary 
objective  of  this  extension  is  to  allow  us  to  acquire  and  process  data  at 
remote  locations  using  interactive  software  similar  to  that  available  in 
our  laboratory. 


We  have  made  considerable  progress  during  the  reporting  interval.  We 
have  defined  a  fonnat--the  microcomputer  flat  file--and  developed  a  number 
of  supporting  programs.  These  include: 


rrsr-. 

1. 

Flat 

file 

keyboard  entry 

2. 

Flat 

file 

digitizer  entry 

3. 

Flat 

file 

editing 

U . 

Flat 

file 

screen  plot 

p. 

5  . 

Flat 

file 

pen  plot 

6. 

Flat 

file 

printer  tabulate 

7  . 

Flat 

file 

select  and  merge 

ft-; 

• 

8. 

Flat 

file 

statistics 

2.4.2  Flat  fiia  transfer  protocols 


In  order  to  analyze  data  acquired  at  a  remote  location  using  a 
microcomputer  such  as  the  APPLE  or  IBM-PC,  it  is  necessary  to  transfer 
floppy  disk  files  to  larger  computers.  During  the  last  year  we  have 
acquired  and  tested  protocols  from  APPLE  lie  to  IBM-PC,  APPLE  He  to 
IBM  3033,  and  IBM-PC  to  IBM  3033.  In  iddition,  we  have  made  some  progress 
in  developing  an  APPLE-IIe  to  HP-i.000  protocol.  Work  is  also  in  progress 
to  allow  transfer  from  the  IBM-PC  to  VAX-'80. 


2.4.3  High  level  software 


Since  our  last  report,  we  have  made  progress  in  transferring  some  of 
our  analysis  programs  from  the  IBM  3033  to  a  VAX  780.  An  important  element 
in  the  transfer  has  been  the  acquisition  for  the  VAX  of  the  SPEAKEZ 
language,  a  very  high-level  interpretative  language  for  management  and 
analysis  of  structured  data.  We  are  currently  working  on  the  creation  of 
additional  modules  to  tailor  the  language  to  our  own  needs.  For  example,  a 
Fast  Fourier  Transform  function  has  being  added.  As  another  example,  we 
have  also  acquired  the  data  management  system  INGRES  and  will  attempt  to 
interface  this  file  manager  to  the  SPEAKEZ  language  in  a  user-transparent 
fashion. 
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